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ABSTRACT 

Wc present ground-based and HST optical observations of the optical transients 
(OTs) of long-duration Gamma Ray Bursts (GRBs) 060729 and 090618, both at a 
redshift of z = 0.54. For GRB 060729, bumps are seen in the optical light curves (LCs), 
and the late-time broadband spectral energy distributions (SEDs) of the OT resemble 
those of local type Ic supernovae (SNe). For GRB 090618, the dense sampling of our 
optical observations has allowed us to detect well-defined bumps in the optical LCs, as 
well as a change in colour, that are indicative of light coming from a core-collapse SN. 
The accompanying SNe for both events are individually compared with SN1998bw, a 
known GRB-supernova, and SN 19941, a typical type Ic supernova without a known 
GRB counterpart, and in both cases the brightness and temporal evolution more 
closely resemble SN1998bw. We also exploit our extensive optical and radio data for 
GRB 090618, as well as the publicly- available SWIFT-XEI data, and discuss the 
properties of the afterglow at early times. In the context of a simple jet-like model, the 
afterglow of GRB 090618 is best explained by the presence of a jet-break at t — to > 0.5 
days. We then compare the rest-frame, peak V-band absolute magnitudes of all of the 
GRB and X-Ray Flash (XRF)-associated SNe with a large sample of local type Ibc 
SNe, concluding that, when host extinction is considered, the peak magnitudes of the 
GRB/XRF-SNe cannot be distinguished from the peak magnitudes of non-GRB/XRF 
SNe. 

Key words: gamma-ray burst: (individual: GRB 060729) — gamma-ray burst: (indi- 
vidual: GRB 090618) — supernovae: general. 



1 INTRODUCTION 

Compelling evidence eoiinceting the generation of bursts of 
gamma-rays and flashes of X-rays with the gravitational 
core-collapse of stripped-envelope supernovae (SNe) contin- 
ues to grow, e.g., XRF 100316D and SN2010bh (Chornock 
et al. 2010; Starhng et al. 2010). The connection between 
GRB 980425 and unusual type Ic SN1998bw (Galama et 
al. 1998) provided an intriguing clue to the so-called "GRB- 
SN" connection (e.g., Woosley & Bloom 2006) that was later 
confirmed with the "smoking-gun" detection of GRB 030329 
and SN2003dh (Stanck et al. 2003; Hjorth et al. 2003; Math- 
eson et al. 2003). 

Supporting the spectroscopic detections of long- 
duration GRBs (i.e., burst durations greater than ^ 2 s) and 
XRF-associated SNe are numerous photometric detections, 
e.g., XRF 020903 (Bersier et al. 2006); GRB 041006 (Stanek 
et al. 2005). Late-time deviations from the power-law de- 
cline expected for a GRB afterglow, that are accompanied 
by a change in colour, are interpreted as evidence of super- 
novae (e.g., Zeh et al. 2004). This spectroscopic and pho- 
tometric evidence appears to strongly favour massive-star 
models, such as the "Collapsar" model (Woosley 1993) and 
the millisecond-magnetar model (e.g., Usov 1992; Thomp- 
son 1994; Wheeler et al. 2000; Zhang & Meszaros 2001), for 
GRB production. 

Indirect evidence supporting the "GRB-SN" connection 
was the realization that GRBs and type Ibc SNe occur in 
the brightest regions of their host galaxies (Fruchter et al. 
2006; Kelly et al. 2008). These observations link GRBs and 
XRFs to sites of star-formation, and consequently to massive 
stars. Similar observations of the positions of Wolf-Rayet 
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(WR) stars, a possible progenitor of GRBs and XRFs, also 
show that their distribution in their host galaxies favour the 
brightest regions (e.g., Leloudas et al. 2010). Additionally, 
observations of GRB host galaxy morphology reveal a high 
fraction of merging and interacting systems (e.g., Conselice 
et al. 2005; Wainwright et al. 2007), which was shown to be 
an efficient trigger of star formation in galaxies (e.g., Joseph 
et al. 1984; Kennicutt & Keel 1984). 

The connection between long-duration GRBs and 

bright supernovae is not without complexity, however. While 
questions surround the nature of GRB 980425 in its relation 
to "classical" GRBs, two nearby events have challenged the 
idea that all long-duration GRBs are accompanied by an 
optically-bright SN. GRBs 060505 and 060614 (Gal- Yam et 
al. 2006; Fynbo et al. 2006; DeUa Valle et al. 2006; Ofek et al. 
2007) had no accompanying supernova down to limits more 
than 100 times fainter than SN1998bw (though the classi- 
fication of GRB 060614 as a long-duration GRB is ques- 
tionable (e.g., Gehrels et al. 2006; Norris et al. 2010). Thus 
the homogeneity of the "GRB-SN" connection, where the 
nai've expectation that all long-duration GRB events pro- 
duce enough nickel to power an optically-bright SN, may be 
an over-simplification of a much more complex reality. 

So while the "GRB-SN" connection has been estab- 
lished, many questions still remain such as: "What kind of 
progenitors produce these events? Are the progenitors all 
the same? Why do some massive stars become GRB/XRFs 
while most do not?" General understanding of the types of 
progenitors that can give rise to a GRB, as well as the en- 
vironments they occur in, have provided indirect clues. The 
occurrence of "optically dark bursts" (i.e., GRBs without 
an optical afterglow) and the (sometimes) large amounts of 
dust in the vicinity of GRB progenitors that can obscure 
light from a supernova is a first attempt to explain the lack 
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of SNe for some nearby GRBs, e.g., GRB 090417B (Holland 
et al. 2010). Additionally, nearby extremely faint (more than 
100 times fainter than SN1998bw) core-collapse SNe have 
been detected (Zampieri et al. 2003; Valenti et al. 200!|3) 
with low expansion velocities and extremely low nickel pro- 
duction. These under-luminous events provide another plau- 
sible explanation (e.g., Tominga et al. 2007; Moriya et al. 
2010) for the lack of accompanying SNe with long-duration 
GRBs as well as illustrating the diverse optical properties of 
core- collapse SNe. 

The evidence gathered so far shows that at least 
some long-duration GRBs and XRFs are accompanied 
by an optically-bright type Ic SN. Individual cases of 
spectroscopically-identified connections of type Ibc super- 
novae to GRBs and XRFs are the ideal way to further our 
understanding, however this is only technologically viable for 
low-redshift bursts. Circumstantial evidence provided pho- 
tometrically has also furthered our understanding by high- 
lighting the enigmatic differences in GRB/XRF-associated 
SNe. 

Here we present results of data obtained on HST and 
ground-based facilities of GRB 060729, and ground-based 
data for GRB 090618, both of which are at a common red- 
shift of 2: = 0.54. In section 2 we present our photometry 
and results for GRB 060729, effectively exploiting our HST 
data to obtain image-subtracted magnitudes that provide 
evidence for an associated supernova. In section 3 we present 
our photometry and results for GRB 090618, which resulted 
in densely-sampled Rc and i band light curves (LCs) that 
clearly show "bumps" that are accompanied by a change 
in colour that we attribute to flux coming from a stripped- 
envelope, core-collapse supernova. In section 4 we compare 
our results of the rest-frame, absolute magnitudes of the SNe 
with those of: (1) aU GRB/XRF-associated SNe and, (2) a 
sample of local type Ibc SNe, concluding that the progeni- 
tors that give rise to GRB/XRF SNe may be similar to those 
that produce local type Ibc SNe but with some differences. 

Throughout the paper, observer-frame times are used 
unless specified otherwise in the text. The respective de- 
cay and energy spectral indices a and /? are defined by 
fuOc{t — to)~"v~^ , where to is the time of burst and v 
is the frequency. Foreground reddening has been corrected 
for using the dust maps of Schlegel et al. (1998), from 
which we find E{B -V) = 0.055 mag for GRB 060729 and 
E[B -V) = 0.085 mag for GRB 090618. We adopt a flat 
ACDM cosmology with a Hubble parameter of Hq — 71 
km/s/Mpc, a matter density of Hm = 0.27, and a cosmo- 
logical constant of JIa = 1 — ~ 0.73. For this cosmology 
a redshift of z — 0.54 corresponds to a luminosity distance 
of dL = 3098.7 Mpc and a distance modulus of 42.45 mag. 



2 GRB 060729 

The Swtft Burst Alert Telescope (BAT) detected the long- 
duration (T90 ~ 115 s; where T90 ) GRB 060729 on July 
29, 2006 (Grupe et al. 2006) with a remarkably bright and 
long-lasting X-ray afterglow (Grupe et al. 2007; Grupe et al. 

1 Though see Foley et al. (2009) and Foley et al. (2010) who 
present various lines of evidence that SN 2008ha was the result 
of a thermonuclear explosion of a carbon-oxygen white dwarf. 



2010). The redshift was measured by Thone et al. (2006) to 
be 2: = 0.54, and a later spectroscopic analysis by Fynbo et 
al. (2009) measured the redshift to be z = 0.5428. 

2.1 Observations & Photometry 

2.1.1 Ground-based Data 

We obtained data with four ground-based telescopes: 
Panchromatic Robotic Optical Monitoring and Polarimetry 
Telescopes (PROMPT), Cerro Tololo Inter- American Obser- 
vatory (CTIO), Faulkes Telescope South (FTS) & Gemini- 
South. Two epochs of data were taken with one of the 0.41m 
Panchromatic Robotic Optical Monitoring and Polarimetry 
Telescopes (PROMPT) ~ 0.5 & 1.5 days after the burst. 
The first epoch yielded images in B and z, while the second 
epoch was only in z. One epoch of data was obtained on the 
Cerro Tololo (CTIO) Blanco 4m telescope 17.6 days after 
the burst, yielding an Rc image. Three epochs of data were 
obtained on the Gemini-South (Gemini-S) 8.1m telescope, 
the first one was 4.6 days post-burst in griz, and two addi- 
tional epochs in r at 46.5 and 49.5 days post-burst. Finally, 
one epoch of images taken was obtained on 5**" April, 2007 
(~ 250 days after the burst) on the 2m Faulkes Telescope 
South (FTS). Observations were made of the GRB field (well 
after the GRB and SN had faded below the instrument de- 
tection limit) in BVRci, as well as images taken of Landolt 
photometric standard regions (Landolt 1992) in the same 
filters. The observations taken by FTS that are used in our 
calibration were made under photometric conditions. 

Aperture photometry wasperformed on all images us- 
ing standard routines in IRAFIj. A small aperture was used, 
and an aperture correction was computed and applied. The 
aperture-corrected, instrumental magnitudes were then cal- 
ibrated via standard star photometry into magnitudes in 
BVRcIc- Using the images of Landolt standards, IRAF rou- 
tines were used to solve transformation equations of the 
form: 

mir,st ^ Zp + M + aiX + a2{B - V) (1) 

where rriinat is the aperture-corrected, instrumental magni- 
tude, M is the standard magnitude, Zp the zero-point, a\ the 
extinction, X the airmass, 02 the colour-term and (B — V^) 
(and variations thereof) the colour. 

The validity of the transformation equations was 
checked by using the solutions on the Landolt standard stars 
in the CCD images, which revealed computed magnitudes 
that were compatible with those in the Landolt Catalogue 
(1992) within the magnitude errorbars (typical a ~ 0.04 
mag). 

The solutions from the standard star photometry were 
then applied to a sequence of secondary standards in the 
field of GRB 060729. The Gemini-S observations are cal- 
ibrated against these stars using a zero-point and colour 
term in filters gRcIcZ. The g and z magnitudes of secondary 
standards in the GRB field were calculated using transfor- 
mation equations from Jordi et al. (2006), and r and i were 

IRAF is distributed by the National Optical Astronomy Obser- 
vatory, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Table 1. Ground-based Photometry of GRB 060729 



T-To (days) 


Telescope" 


Filter 


Mag 


a (mag) 


Calibrated to 


0.452 


JT 1 Olll L) L 




18.22 


0.09 


Vega 




Prompt 






n OR 
u.uo 


AB 


0.481 


Prompt 


z 


17.29 


0.07 


AB 


0.515 


Prompt 


z 


17.22 


0.06 


AB 


0.551 


Prompt 


z 


17.49 


0.07 


AB 


0.571 


Prompt 


B 


18.40 


0.04 


Vega 


0.594 


Prompt 


z 


17.49 


0.07 


AB 


1.574 


Prompt 


z 


18.43 


0.17 


AB 


4.601 


Gemini-S 


Rc 


20.24 


0.04 


Vega 


4.613 


Gemini-S 


a 


20.91 


0.04 


AB 


4.626 


Gemini-S 


Ic 


19.92 


0.03 


Vega 


4.638 


Gemini-S 


z 


20.51 


0.04 


AB 


17.598 


CTIO 


Rc 


21.90 


0.06 


Vega 


46.583 


Gemini-S 


Rc 


23.29 


0.06 


Vega 


49.554 


Gemini-S 


Rc 


23.32 


0.05 


Vega 



"Telescope key: CTIO: 4m Cerro Tololo Blanco Telescope; Prompt: 0.41m Prompt Telescope; Gemini-S: 8.1m Gemini-South Telescope. 

All magnitudes have been corrected for foreground extinction. 



Table 2. HST Photometry of GRB 060729 



T - To (days) 


Filter" 


Mag'' 


cr (mag) 


Calibrated 


8.983 


F330W 


22.00 


0.03 


AB 


17.051 


F330VK 


23.40 


0.05 


AB 


26.125 


F330VK 


24.08 


0.10 


AB 


9.042 


F625VK 


21.09 


0.01 


Rc (Vega 


16.909 


F625W 


21.94 


0.02 


Rc (Vega 


26.250 


F625W 


22.57 


0.03 


Rc (Vega 


48.689 


F625W 


24.07 


0.07 


Rc (Vega 


139.888 


F625VF 


26.20 


0.62 


Rc (Vega 


9.059 


F850LP 


20.90 


0.01 


Ic (Vega) 


16.926 


F850LP 


21.74 


0.02 


Ic (Vega) 


26.266 


F850LP 


22.24 


0.03 


Ic (Vega) 


48.757 


F850LP 


23.28 


0.05 


Ic (Vega) 


139.905 


P850LP 


25.65 


0.57 


Ic (Vega) 


9.139 


F160W 


21.14 


0.02 


AB 


16.500 


F160W 


22.16 


0.04 


AB 


27.198 


F160W 


22.85 


0.05 


AB 


9.124 


F222M 


20.38 


0.30 


AB 


16.460 


F222M 


> 22.1 




AB 


27.127 


F222M 


> 21.7 




AB 



"Filter key: F330W, F625W, F850LP: ACS; FWOW, F222M: 
NICMOS. 

''Apparent magnitudes are of the OT from the subtracted images 
except for the ACS F330W and NICMOS F222M filters (for 
which no image subtraction occurred). 

All magnitudes have been corrected for foreground extinction. 



2.1.2 Hubble Space Telescope Data 

We obtained a total of six epochs of data with HST, us- 
ing the Advanced Camera for Surveys (ACS) and the High- 
Resolution Channel (HRC) in F330W; and the Wide Field 
Channel (WFC) in F625W and F850LP; the Near In- 
frared Camera and Multi-Object Spectrograph (NICMOS) 
in F160W and F222M; and the Wide-Field Planetary Cam- 
era 2 (WFPC2) in F300W, F622W and i^850LP for the last 
epoch. The images were reduced and drizzled in a standard 
manner using Multidrizzle (Koekemoer et al. 2002). The use 
of WPFC2 images in the last epoch was a consequence of 
the ACS being broken at this time. While we were wary 
that images from two different cameras and two different 
filters might disrupt our intention of using the last epoch 
as a template image for use in image subtraction, we were 
able to work around the problem of the differing cameras 
and filters by carefully choosing the stars used in the image- 
subtraction procedure (see section 2.1.3). The ACS images 
that were used for image subtraction were drizzled to the 
WFPC2 scale (0.10"), while the remaining images were driz- 
zled to the native scale of the individual cameras (0.2" and 
0.10" for NICMOS and WFPC2 respectively). 

Image subtraction was not carried out on the F330W 
images because the host is not visible in the final WFPC2 
{F300W) image, thus we used the procedure described in 
Appendix B of Sirianni et al. (2005) for the ACS data: (a) 
aperture photometry using aperture of 0.15", (b) aperture 
correction for apertures from 0.15" to 0.50", (c) aperture 
correction for F330W for aperture 0.50" to infinity (Table 
5 of Sirianni et al. 2005), (d) CCD charge transfer efficiency 
(CTE), (e) used AB zeropoint. 



calibrated to Rc and Ic. The single CTIO epoch was cal- 
ibrated directly against the stars in the GRB field using 
a zero-point (no colour term). Images taken by PROMPT 
are calibrated against standard stars observed by PROMPT 
in B and z, using a zero-point and colour term. All of the 
ground-based photometry, which has been corrected for fore- 
ground extinction, is presented in Table [1] 



2.1.3 Image Subtraction 

Using ISIS (Alard 2000) we performed image subtraction on 
the ACS F625W and _F850LP images, using the WFPC2 
images in F622W and F850LP as the respective templates. 
The ACS images were drizzled to the WFPC2 scale (0.10") 
prior to subtraction. Image subtraction was also performed 
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Table 3. HST Photometry of the Host Galaxy of GRB 060729 



T - To (days) 


Filter" 


Mag 


a (mag) 


Calibrated to 


425.298 


F300VK 


> 24.80 




AB 


428.984 


F622W 


24.42 


0.20 


Rc (Vega) 


429.181 


F850LP 


24.26 


0.20 


Ic (Vega) 


425.363 


Fieow 


22.82 


0.06 


AB 


424.806 


F222M 


> 22.10 




AB 



"Filter key: F622Vl^ & F850LP: WFPC2; F160W: NICMOS. 
All magnitudes have been corrected for foreground extinction. 



on the NICMOS F160W images, using the last epoch as the 
template. Subtracting the WFPC2 epoch from the first five 
ACS epochs gave a clear detection of the optical transient 
(OT). As the reference image was taken with a different 
camera and different filters (with the corresponding differ- 
ences in the filters' transmission curves and different CCD 
responses) we checked for the presence of a colour-term in 
the residuals in the subtracted images. For very red objects 
{F622W - F850LP > 1.5) a smaU colour dependent effect 
was seen, however when we restricted the colour range to in- 
clude only stars with < F622W - F850LP sC 1.0 we found 
no statistically-significant colour-term in the residuals. The 
colour of the OT was never beyond this colour range, thus 
by using the stars in this colour range to match the ACS 
frames to the WFPC2 template in our image-subtraction 
procedure we are confident of the validity of our results. 

Aperture photometry was then performed on the sub- 
tracted images using an aperture of 0.15", and correct- 
ing for steps (a) - (c) as described previously but for the 
WFPC2 images (Holtzman et al. 1995). The instrumental 
magnitudes were then calibrated to Rc and 7c using a zero- 
point and colour-term. Aperture photometry was performed 
on the NICMOS F160W images using the procedures pre- 
scribed in the NICMOS data handbook (Thatte et al. 2009). 
We used an aperture of 0.40" and corrected for: (a) aperture 
correction for apertures from 0.40" to 1.00", (b) correct to 
infinite aperture by multiplying the flux (counts) by 1.075, 
(c) used AB zeropoint. All of the HST magnitudes, which 
have been corrected for foreground extinction, are listed in 
Table m 



2.1.4 Host Photometry 

The host was clearly visible and extended in all of the fl- 
nal epoch images taken by WFPC2, apart from the F300W 
image. Aperture photometry was performed on the images: 
WFPC2 {F622W and F850LP) and NICMOS (F160W), 
using an aperture of 0.7", 0.4" and 0.6" respectively. As the 
host is clearly extended, CTE effects have been neglected 
when determining the host magnitudes as the effects due to 
CTE on WFPC2 images is only well understood for point- 
sources. The corrected instrumental magnitudes in F622W 
and F850LP were then calibrated to Rc and Jc using a zero- 
point and colour-term, while the F160W and F222M mag- 
nitudes are given as AB magnitudes. All of the host magni- 
tudes, which have been corrected for foreground extinction, 
are listed in Table [S] 



2.2 Results &: Discussion 

2.2.1 The Afterglow 

Grupe et al. (2007) found that the optical and X-ray light 
curves of GRB 060729 displayed similar temporal behaviour 
(i.e., the decay constants and time of break in the X-ray, UV 
and optical filters were consistent up to late times; see Figure 
8 in Grupe et al. 2007). Thus, assuming at late times that the 
decay constant will be approximately the same in all optical 
passbands, we used our HST{ACS) F330W magnitudes to 
characterise the temporal behaviour of the GRB afterglow 
at longer wavelengths (i.e., Rc and 7c). The reason for using 
the F330W magnitudes is that any SN contribution at these 
wavelengths is expected to be negligible with respect to the 
afterglow (e.g., see Figure [T}. 

Using the Ultraviolet and Optical Telescope (UVOT) 
aboard Swift, Grupe et al. (2007) found that the optical 
afterglow was visible in the UVOT-?7 filter up to 7.66 days 
after the burst. We clearly detect the GRB afterglow up 
to our third HST epoch (26.12 days after the burst) in the 
HST{ACS) F330W fiher. Thus, by combining our HST data 
with that obtained by Swift, we were able to extend the time- 
domain of the analysis and more accurately determine the 
temporal behaviour of the afterglow. 

To enable us to compare the F330W magnitudes ob- 
tained with HST (ACS) to those determined by Swift, we 
used IRAF/Synphot/Calcphot to transform our 7^33014^ 
magnitudes into U and corrected for foreground extinction. 
Next, using transformation equations determined by Swift 
(Poole et al. 2008), we transformed the UVOT-C/ magni- 
tudes into U. All of the U-hand data {HST and Swift) were 
thus calibrated to U. 

A broken power- law was fitted to the data (Figure [TJ, 
using the form: 



m„{t) 



-2.5x 



+ 



Th, 



-73(2) 



where t is the time since the burst and Threak is the time 
when the power- law changes from temporal index qi to 02. 

The values of the fit <od jdof = 1.31) are: ai = 0.01 ± 
0.03; Q2 = 1.65 ± 0.05 and Tbrcak = 0.75 ± 0.08 days. When 
we restrict the upper limit to the time range to 8 days (i.e., 
up to the last UVOT-f/ detection) we flnd {x^/dof = 1.33): 
02 = 1.47 ± 0.11 and Tbrcak = 0.58 ± 0.10 days, which are 
consistent with the decay constant and time of break found 
by Grupe et al. (2007) for the UVOT-[/ data (see Table 5 
of Grupe et al. 2007). We note that we originally fit the 
data with the complete Beuermann Function (Beuermann 
et al. 1999) and let n vary (where n is a measure of the 
"smoothness" of the transition from ai to 02). The results 
of our fit gave n « 1 (i.e., a smooth transition), so we fixed 
n — 1 when performing the final fit. 

In Figure[l] variability is seen in the early-time data. To 
check that the early-time data does not affect our empirical 
fit at late-times, we fit a single power-law to the data at 
t~to > 2.0 days. It was seen that by limiting the fit to data 
after this time, we found a = 1.64 ± 0.07, fully consistent 
with the value of 02 found previously. 

The value of 02 = 1.65 ± 0.05 is consistent with that 
found by Grupe et al. (2010) for the X-ray light curve 



(a 



1.61 



-1-0.10 

0.06 



) from 1.2 X 10" s (13.8 days) to 38 x 10" 



s (439.8 days), though we note that the optical and X-ray 
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Figure 1. GRB 060729: Light-curve of UVOT (U) (triangles) 
and HST (F330W) (squares) magnitudes transformed into U. 
The afterglow model is shown as a broken power-law (solid line), 
where oi = 0.01±0.03; 02 = 1.65±0.05 and Tt^eafe = 0.75±0.08 
days. Plotted for comparison (pentagons) is the UVOT (UVWl) 
LC of SN 2006aj/GRB 060218 (Brown et al. 2009) as it would 
appear at z = 0.54. All magnitudes have been corrected for fore- 
ground extinction and the LC of 2006aj has also been corrected 
for host extinction. At these wavelengths the contribution of flux 
from an accompanying SN is expected to be negligible in compar- 
ison to the afterglow. This assumption is confirmed by the HST 
observations. 
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Figure 2. GRB 060729: afterglow-subtracted, supernova Rc 
(blue) and Ic (black) light curves. The solid points are derived 
from the image-subtracted magnitudes from the HST data while 
the open points are derived from the magnitudes obtained by 
mathematically subtracting the host flux from the Gemini-S and 
CTIO epochs (see text). Plotted for comparison are the LCs 
of SN1998bw (dashed) and SN1994I (dot-dashed) in the same 
colours, as they would appear at 2 = 0.54. The brightness and 
time-evolution of the SN resembles that of SN1998bw rather than 
SN1994I. All magnitudes have been corrected for foreground and 
host extinction. 



decay constants do not necessarily have to be the same, as 
seen for many bursts (e.g., Melandri et al. 2008). 

Plotted for comparison in Figure [T] is the UVOT 
(UVWl) LC of XRF 060218/SN 2006aj. At a redshift of 
z = 0.54, the light emitted in the rest-frame at UVOT 
(UVWl) wavelengths will be redshifted into observer- frame 
U. We can see that the broken power-law fits the U data well 
and we conclude that the contribution of flux at this wave- 
length from an accompanying SN is negligible in comparison 
to the flux from the afterglow. 



2.2.2 The Supernova 

For the complete GRB event, we assume that flux is com- 
ing from up to three sources: the afterglow, the host and a 
(possible) SN. We have already photometrically removed the 
flux due to the host in the HST images via image subtrac- 
tion, however we were not able to perform image subtraction 
on the Gemini-S and CTIO images. Instead we mathemati- 
cally subtracted the contribution of host flux from the three 
Gemini-S and CTIO images, after all of the magnitudes were 
converted into fluxes using the zero-points from Fukugita et 
al. (1995). Thus, after removing the flux due to the host, and 
then subtracting the flux due to the GRB afterglow (which 
was characterised by Equation 2), any remaining flux could 
be due to an accompanying SN. 

Our host- and afterglow-subtracted light curves are 
shown in Figure (2] where the quoted errors have been cal- 
culated in quadrature. Plotted for comparison are the dis- 
tance and foreground-corrected light curves of SN1998bw 
and SN1994I. The SN1998bw LCs are model light curves 
(i.e., interpolated/extrapolated through frequency) however 



we have used the rest-frame B and V LCs of SN1994I to 
represent the observer-frame Rc and Ic LCs respectively. It 
is seen, even despite the paucity of data-points, that the SN 
associated with GRB 060729 resembles more closely that of 
SN1998bw than SN1994I in regards to both time evolution 
and brightness. 

We find peak apparent magnitudes of the associated SN 
of Rc = 23.80 ± 0.08 and Ic = 23.20 ± 0.06. To make an ac- 
curate estimate of the peak absolute magnitude of the SN, 
we have determined the amount of host extinction (Section 
2.2.4) to be E{B — V) — 0.61 mag. In comparison, Grupe 
et al. (2007) determined from their optical to X-ray spec- 
tral energy distribution (SED) analysis of the afterglow a 
rest-frame extinction local to the GRB of E{B — V)= 0.34 
mag. However, using the same UVOT data-set as Grupe et 
al. (2007) as well as additional Rc data, Schady et al. (2010) 
found a rest-frame extinction of E{B — V) ^ 0.06. We note 
that our value of the host reddening is poorly constrained by 
our SED template fitting (though it is worth noting that the 
best two fits to the host SED both imply large amounts of 
reddening). However, it has been seen that the host extinc- 
tion for many bursts is not always the same as the extinction 
local to the event (e.g., GRB 000210 (Gorosabel et al. 2003), 
where the local reddening was very high, but the host ex- 
tinction was negligible. Further examples are seen in Modjaz 
et al. 2008 and Levesque et al. 2010). Thus it is more rea- 
sonable to use the extinction derived by Schady et al. (2010) 
as it is local to the event and is a more complete analysis 
than that performed by Grupe et al. (2007). 

As the majority of long-duration GRB host galaxies 
contain dust that is Small Magellanic Cloud (SMC)-like 
(e.g., StarUng et al. 2007; Schady et al. 2007; Kann et al. 
2010), we applied the SMC extinction law (Ry = 2.93; Pei 
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Figure 3. GRB 060729: Observer-Frame SED. The host- 
subtracted HST data (filled circles) of the OT (afterglow and 
supernova) has been corrected for foreground and host extinc- 
tion. Plotted for comparison are the (see text) UBVRI data of 
SN1998bw at 25.8 days and the UBVRI data for SN1994I at 28.8 
days, which have been corrected for foreground and host extinc- 
tion. At 9.0 days the SED is fitted by a power-law with spectral 
index /3 = 0.59 ± 0.02, however at 26.2 days the SED resembles 
that of the local Ic SNe, although with less curvature. 



1992) to the value oiE{B- V) found by Schady at al. (2010) 
to estimate the rest-frame extinction. We then approximated 
the rest-frame V-band magnitudes from the observer-frame 
7c-band magnitudes, and find an absolute magnitude of the 
SN associated with GRB 060729 of My = -19.43 ± 0.06. 
This implies that the SN associated with GRB 060729 has 
approximately the same peak brightness as 1998bw in the 
V-band (see Table [T]). 



2.2.3 The Spectral Energy Distribution 

We determined the observer-frame spectral energy distri- 
bution of the OT (afterglow and SN) for two epochs: 9.0 
and 26.2 days after the burst (Figure |3]). We used the zero- 
points from Fukugita et al. (1995) to convert the magni- 
tudes to fluxes, and corrected for foreground and host extinc- 
tion. Plotted for comparison are the observer-frame SEDs for 
SN1998bw (25.8 days) and SN1994I (28.8 days), which have 
been corrected for distance as well as foreground and host 
extinction. At 9.0 days the SED is fitted by a power-law with 
spectral index /3 = 0.59 ± 0.02. However, the most striking 
feature of Figure|3]is at 26.2 days, where the SED resembles 
that of 1998bw and 19941 at a similar epoch, although with 
less curvature. 



F850LP = 24.60 ±0.20, F160W^ = 24.12 ±0.0^, and upper 
limits: F330 > 24.80, F222M > 22.10 mag. 

We used the procedure described in Svensson et al. 
(2010) to calculate the stellar masses, star-formation rate 
and the specific star-formation rate of the host galaxy. Us- 
ing a redshift of 2 = 0.54, we fitted our HST magnitudes 
to synthetic galaxy evolution models (Bruzual & Chariot 
2003). None of the synthetic models were very well con- 
strained, however the best fit template (x^/do/ — 6.58) is 
for a galaxy with a young stellar population and low metal- 
licity. We note that the next best-fit model is for a M82-type 
galaxy with large amounts of reddening. 

The magnitudes of the best-fitting template are: Mu ~ 
-16.72, Mb = -17.13, My = -17.49, Mk = -19.94 
and E{B — V) ~ 0.61 mag. The magnitudes are not 
host-extinction corrected. We also find a star-formation 
rate of log{SFR[MQ/yr]) = -0.89+o;o7; a stellar mass of 
log(Mass[MQ]) = 9.13lQQg and a specific star-formation 
rate of log{S S F R[Gyr-^]) = -1.03l^;o^. The latter values 
of the star-formation rates and mass are estimated from an 
extinction-corrected SED. The quoted errors are la and are 
systematic only and are calculated from the distributions of 
SFR, mass, etc., given by all templates (i.e., the errors re- 
fiect how good the best-fit parameters are compared with 
all of the considered models). 

We were able to determine the position of GRB 060729 
in the host galaxy using the HST images taken by WFPC2. 
We find for the position of the OT in our subtracted im- 
ages (error = 0.02"): 06''2T"31.77" , -62'^22'" 12.21"; which 
is offset from the apparent centre of the host by 0.33" ±0.02" 
(2.1 ± 0.1 kpc). 

Secondly, we used the pixel statistic created by Fruchter 
et al. (2006), which calculates the fraction of light (Fught) 
contained in regions of lower surface brightness than the re- 
gion containing a GRB or SN (where a value of Fught ~ 
100.0 corresponds to the event occuring in the brightest re- 
gion of the host galaxy). For GRB 060729, Fught = 39.0 in 
the FQ22W filter, which is somewhat smaller than that seen 
in other GRBs by Fruchter et al. (2006). 

Thus we conclude that the host galaxy of GRB 060729 
has a young stellar population with a modest stellar mass (~ 
10^ M0) and star-formation rate {SFR ~ 0.13 MQyr"^), 
suggesting that the host is typical of other GRB host galax- 
ies (e.g., Wainwright et al. 2005; Fruchter et al. 2006; 
Savaglio et al. 2009; Svensson et al. 2010; Christensen et 
al. 2004). 



3 GRB 090618 

The Swift Burst Alert Telescope (BAT) detected the long- 
duration (T90 = 113. 2±6 s) GRB 090618 on June 18, 2009 at 
08:28:29 UT (Schady et al. 2009) which was followed up by 



2.2.4 The Host Galaxy 

The detection of the host galaxy with WFPC2 allowed us to 
construct its optical SED (Figure^. Magnitudes of the host, 
which are corrected for foreground reddening, were trans- 
formed to AB magnitudes, giving: _F622W^ = 24.54 ± 0.20, 



^ A recent image of the host galaxy has been taken with the Wide 
Field Camera 3 aboard HST in filter F160VK by A. Levan, who 
confirms that the host brightness in this filter has not changed 
between this epoch and our epoch taken at 425 days (A. Levan, 
private communication). 



8 Z. Cano et al. 



10' r- 




Figure 4. GRB 060729: SED of the host galaxy. The best-fit 
model is for a galaxy with a young stellar population and low 
metallicity. 

many teams with ground-based telescopes at optical wave- 
lengths, and published in the GCN CircularsQ 

The redshift was measured to be z = 0.54 by Cenko et 
al. (2009) and FatkhuUin et al. (2009). 

3.1 Observations & Photometry 

3.1.1 Optical Data 

We obtained optical observations with many ground-based 
telescopes: the 2m Faulkes Telescope North; the 2m Liver- 
pool Telescope; the 8.1m Gemini-North; the 4.2m William- 
Herschel Telescope (WHT); the 2.5m Isaac Newton Tele- 
scope (INT); the 2.6m Shajn telescope of CrAO; the 1.25m 
AZT-11 telescope of CrAO; the 1.5m Maidanak AZT-22 
telescope with the SNUCAM CCD camera (Im et al. 2010); 
the 1.5m AZT-33IK telescope of Sayan observatory, Mondy; 
the 6m BTA-6 optical telescope and the Im Zeiss- 1000 tele- 
scope of the Special Astrophysical Observatory of the Rus- 
sian Academy of Sciences (SAO-RAS); the Im LOAO tele- 
scope (Lee, Im & Urata 2010) at Mt. Lemmon (Arizona, 
USA) which is operated by the Korea Astronomy Space Sci- 
ence Institute; 15.5cm telescope at Bohyunsan Optical As- 
tronomy Observatory (BOAO) in Korea which is operated 

^ The data presented in this paper supersedes those in the GCN 
Circulars: FTN: Melandri et al. 2009 and Cano et al. 2009; LOAO: 
Im et al. 2009a; OAUV-0.4m telescope: Fernandez-Soto et al. 
2009; Shajn: Rumyantsev & Pozanenko 2009; BOAO: Im et al. 
2009b; SAO-RAS: FatkhuUin et al. 2009 (photometry & spectra); 
HOT: Anupama, Gurugubelli & Sahu 2009; Mondy: Klunko, Vol- 
nova & Pozanenko 2009. GRB 090618 was also detected at ra- 
dio wavelengths and published in GCN Circulars: AMI: Pooley 
2009a, 2009b; VLA: Chandra & Frail 2009; WSRT: Kamble, van 
der Horst & Wijiers 2009. 



Table 4. Ground-based Optical Photometry of GRB 090618" 



T- To (days) Filter Mag*" cr (mag) Telescope'^ 



0.057 


B 


17.34 


0.05 


FTN 


0.059 


B 


17.38 


0.05 


FTN 


0.062 


B 


17.43 


0.05 


FTN 


0.067 


B 


17.50 


0.06 


FTN 



"Only a section of the data is shown here. The complete table is 
available online at: 

''The apparent magnitude of the GRB-f SN-fHOST. 
'^Telescope key: AZT-11: 1.25m AZT-11 Telescope of CrAO; 
BOAO: the 15.5cm BOAO Telescope; BTA-6: 6m SAO-RAS Tele- 
scope FTN: 2m Faulkes Telescope North; Gemini (N): 8.1m Gem- 
ini North; HCT: 2m Himilayan Telescope; INT: 2.5m Isaac New- 
ton Telescope; LOAO: 1.0m LOAO Telescope; LT: 2m Liverpool 
Telescope; Maidanak: 1.5m Maidanak Telescope; Mondy: 1.5m 
AZT-33IK Telescope; OAUV: 0.4m OAUV Telescope; Shajn: 
2.6m Shajn Telescope of CrAO; WHT: 4.2m William Herschel 
Telescope. 

by the Korea Astronomy Space Science Institute; the 2m Hi- 
malayan Chandra Telescope (HCT) and the T40 telescope 
at the Observatorio de Aras de los Olmos, operated by the 
Observatori Astronomic de la Universitat de Valencia. 

Aperture photometry was performed on all images us- 
ing standard routines in IRAF. A small aperture was used, 
and an aperture correction was computed and applied. The 
aperture-corrected, instrumental magnitudes of stars in the 
GRB field were calibrated to the SDSS catalogue (in filters 
BVRci, where non-SDSS magnitudes were calculated us- 
ing transformation equations from Jordi et al. (2006) using 
Equation 1, initially using a zero-point and a colour term. 
However, it was found that the contribution of a colour term 
to the calibration was negligible, hence only a zero-point was 
used in all cases. All optical magnitudes have been corrected 
for foreground extinction and are listed in Table H) 

3.1.2 X-ray Data 

The X-ray afterglow was detected and monitored by Swift- 
XRT (Schady et al. 2009; Evans et al. 2009; Beardmore & 
Schady 2009). The XRT data have been processed with the 
HEASOFT package v. 6. 9 and the corresponding calibration 
files (standard filtering and screening criteria have been ap- 
plied). For t < 240 s the data are strongly affected by pile- 
up: piled-up Window Timing (WT) mode data have been 
corrected by eliminating a strip of data from the original 
rectangular region of event extraction: the strip size has 
been estimated from the grade distribution (according 
to the Swift nomenclature: see Burrows et al. 2005; Ro- 
mano et al. 2006) . Piled-up Photon Counting (PC) data have 
been extracted from an annular region whose inner radius 
has been derived comparing the observed to the nominal 
point spread function (PSF, Moretti et al. 2005; Vaughan 
et al. 2006). The background is estimated from a source- 
free portion of the sky and then subtracted. The 0.3-10 keV 
background-subtracted, PSF- and vignetting-corrected light 
curve has been re-binned so as to ensure a minimum num- 
ber of 25 re-constructed counts. The count-rate light curves 
are calibrated into flux and luminosity light curves using a 
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time dependent count-to-flux conversion factor as detailed 
in Margutti et al. (2010). In this way the strong spectral 
evolution detected at t < 240 s is properly accounted for. 



3.1.3 Radio Data 

Radio observations were performed with three radio tele- 
scopes: the Arcminute Microkelvin Imager (AMI) Large Ar- 
ray, part of the Cavendish Astrophysics Group, MuUard Ra- 
dio Astronomy Observatory, Cambridge UK; the 22m RT-22 
radio telescope, in Ukraine; and the Westerbork Synthesis 
Radio Telescope (WSRT), in the Netherlands. 

One epoch was obtained on the AMI on 19 June 2009 in 
the 14.6 to 17.5 GHz band at 0.823 days post burst which led 
to a positive detection, reported in the GCN Circulars (Poo- 
ley 2009a). The typical noise level of the maps was 0.115 to 
0.190 mjy, and resolution was typically 40 x 20 arcseconds. 

Three epochs were obtained on the RT-22 radio tele- 
scope: 18 June 2009 (22 GHz), 21 June 2009 (2 GHz and 
8 GHz) and 28 June 2009 (2 GHz and 8 GHz). The GRB 
was not detected in any observation, leading to upper limits 
only. 

Additional radio observations were made with WSRT 
at 2.3 GHz and 4.8 GHz. We used the Multi Frequency 
Front Ends (Tan 1991) in combination with the IVC-I-DZB 
back encjfl in continuum mode, with a bandwidth of 8x20 
MHz. Gain and phase calibrations were performed with the 
calibrators 3C 48 and 3C 286, at 2.3 and 4.8 GHz, respec- 
tively. The observations have been analysed using the Mul- 
tichannel Image Reconstruction Image Analysis and Display 
(MIRIAD; Sault et al. 1995) software package. 

The first observation at 4.8 GHz, 1.7 days after the 
burst, was reported in Kamble et al. (2009), but a careful 
reanalysis of the data resulted in a lower flux value and lower 
uncertainty in that value. All of the radio observations are 
hsted in Table [S] 



3.2 Results & Discussion 

3.2.1 The Optical, X-ray & Radio Afterglow 

Over 5 hours of optical data were obtained on FTN in 
BVRci starting 1.3 hours after the initial GRB trigger. This 
was supplemented by data collected on the other telescopes, 
resulting in well-sampled light curves. Figure [S] shows the 
flrst six days of optical data: this includes all of the detec- 
tions in B and V, and the first six days of data in Rc and 
i. Several days after the burst the onset of light from an ac- 
companying SN starts to become significant, thus the data 
after this time cannot be reliably used to characterise just 
the afterglow as it would incorrectly imply a slower afterglow 
decay. 

Inspection of the panchromatic light curve clearly shows 
an achromatic break. Broken power-laws have been simul- 
taneously fit to the BVRci data from 0.03 days to 6 
days using Equation (2)0, resulting in the following values 
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Figure 5. GRB 090618: The X-ray and optical light curve. 
A broken power-law has been fitted separately to the X-ray 
{x^/dof = 1.17: ai = 0.79 ± 0.01; 02 = 1.74 ± 0.04 and 
Tbrcak = 0.48 ± 0.08 days) and optical data (x'^/dof = 1.34: 
ai = 0.65 ± 0.07; 02 = 1.57 ± 0.07 and Ttreak = 0.50 ± 0.11 
days). The dashed grey lines represent the times of the SEDs. 



ix /dof = 1.38): ai = 0.65 ± 0.08; a2 = 1.57 ± 0.07 and 
Ttreak = 0.50 ± 0.11 days. 

The 0.3-10 keV light curve, also shown on Figure (5] 
has been corrected for foreground and host Nh absorp- 
tion [Nh = 4.0 ± 0.9 X 10^^ cm~^, see below) and displays 
the common "steep-shallow-steep" temporal behaviour (e.g., 
Nousek et al. 2006; Zhang et al. 2006; Evans et al. 2009) seen 
for many other GRBs. The LC is quite featureless (i.e., no 
flares, etc), which is typical of other classical GRBs with 
confirmed SNe (e.g.. Figure 7 in Starling et al. 2010). The 
data at times greater than t — to — 0.006 days were fit with 
a broken power-law (Equation 2) resulting in the following 
values ix^/dof = 1.17): ai = 0.79 ± 0.01; 02 = 1.74 ± 0.04 
and Tbreak = 0.48±0.08 days. While the values of 02 derived 
from the optical and X-ray data are relatively similar, they 
differ by ~ 2cj, which could be explained by the presence 
of a frequency break (i.e., cooling break) between the two 
wavelength regimes (see section 3.2.3 and Figure [TT|) . 

We note that the post-break decay index of a < 2 is 
not "permitted" in the most simple theories. However, shal- 
low post-break decays are not a very rare phenomenon, and 
several examples are seen in the literature (e.g. Zeh et al. 
2006). 

The column density along the line of sight to the GRB 
within the host galaxy was estimated from the absorption 
in the late time X-ray spectrum: the Galactic contribu- 
tion, Nh = 5.8 X 10^" cm~^, was fixed to the value mea- 
sured from 21cm line radio surveys (Kalberla et al. 2005), 
while the intrinsic (GRB rest-frame) Nh was found to be 
Nh = 2.7 ± 0.3 x 10^^ cm"^ The Galactic and intrinsic 
Nh were obtained with the tbabs and ztbabs models under 
xspec respectively. Using Nh = 2.7 ± 0.3 x 10^^ cm~^ as 



the rest-frame column density we estimate the amount of 
extinction expected at optical wavelengths using Equation 

^ See sect. 5.2 at |http: / / www. astron. nl / radio-observatory/ astronomers / wsrt-guide-observations | 
^ Similar to our analysis of the afterglow of GRI3 060729, we orig- 
inally fit the optical and X-ray data with the complete Beuermann suits of our fit found n 1, so we fixed n = 1 when performing 
Function (Beuermann et al. 1999), letting n vary. Again, the re- the final fit. 
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Table 5. Radio Observations" of GRB 090618 



T-To (days) 


F2.3 GHz 


O" -P4.8 GHz 


a 


-^8.46 GHz O" -P15 GHz -^22 GHz 


<j Telescope'' 


0.168 








<5 


RT-22 


0.824 








0.620 0.115 


AMI 


0.980 








0.383 0.047 


VLA'^ 


1.680 




0.388 


0.030 




WSRT 


3.203 


<3 






<3 - - - - 


RT-22 


4.670 




0.093 


0.045 




WSRT 


4.700 


0.042 


0.052 






WSRT 


8.680 




0.078 


0.030 




WSRT 


10.258 


<3 






<8.900 _ _ - _ 


RT-22 
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''Telescope key; AMI: AMI Large Array; RT-22: 22-m RT-22 dish; WSRT: Westerbork Synthesis Radio Telescope; VLA: Very Large 
Array. 

Observations from GCN 9533 (Chandra & Frail 2009). 
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Figure 6. GRB 090618: Radio light curves at 2.3 GHz, 4.8 GHz, 
8.46 GHz and 16 GHz, where upper limits are denoted by filled 
triangles. Also plotted are the derived radio light curves for a 
spherical (f„{t) oc i^/^) (dotted) and jet-Uke (/^(t) oc t'^^'^) 
(solid) evolution of the afterglow. 

(4) and Table 2 from Pei (1992) for SMC interstellar dust 
parameters, finding a modest value of Ay ~ 0.29 mag. Sim- 
ilarly low values are also derived for Milky Way (MW) and 
Large Magellanic Cloud (LMC)-type templates. 

The radio detections and upper limits are plotted in Fig- 
ure [G] Also plotted are our predicted radio LCs (see section 
3.2.3). The radio detections are consistent with a jet-like 
evolution of the afterglow into an ISM environment. The 
modelling of the radio detections will be discussed further 
in section 3.2.4. 

3.2.2 The Supernova 

We measured the host magnitudes from our late-time WHT 
images (354.7 days) to be Rc = 23.44 ± 0.06 and i = 
23.22±0.06. The fiux contribution from the host was mathe- 
matically subtracted from all epochs, after all of the magni- 
tudes were converted into fluxes using the zero-points from 
Fukugita et al. (1995). The resulting hght curves, which 
have been corrected for foreground extinction, are displayed 
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Figure 7. GRB 090618: Rc (red) and i (blue) light curves of the 
host-subtracted (see text) magnitudes, which have been corrected 
for foreground extinction. The solid line is our model (afterglow 
and SN1998bw-type SN). The afterglow is modelled with a broken 
power-law with parameters: ai = 0.65 ± 0.07; 02 = 1-57 ± 0.07 
and Ttreak = 0.50 ±0.11 days, and the template SN (SN1998bw) 
is dimmed by 0.75 and 0.5 mag in Rc and i respectively. 



in Figure [71 and the quoted errors have been calculated in 
quadrature. Bumps are clearly detected in both hght curves. 

We then subtracted the fiux due to the GRB from the 
"host-subtracted" Ught curves using Equation (2), produc- 
ing the "SN" light curves shown in Figure [8] Also plotted 
are the distance and extinction (foreground) corrected light 
curves of 1998bw and 19941. As was seen for the SN associ- 
ated with 060729, the time evolution and brightness of the 
SN associated with 090618 resembles more closely that of 
SN1998bw than SN1994I. 

We find for the SN associated with GRB 090618 peak 
apparent magnitudes of Rc = 23.45 ± 0.14 and i = 23.00 ± 
0.09. These values, as well as the peak time, are in good 
agreement with the values predicted by Dado & Dar (2010). 
Using the observer-frame i-band magnitudes to approximate 
the rest-frame l/-band magnitudes, as well as our estimation 
of the rest-frame extinction {Ay ^ 0.3 ±0.1; see section 
3.2.3), we find = -19.75 ± 0.13, which is ~ 0.3 mag 
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Figure 8. GRB 090618: Rc (red) and i (blue) supernova light 
curves of the host and GRB-subtracted (see text) magnitudes. 
Plotted for comparison are the light curves of SN1998bw (dashed) 
and SN1994I (dot-dashed) plotted in the same colours, as they 
would appear at z = 0.54. As seen for the 060729-SN, the 
brightness and time-evolution of the 090618-SN resembles that 
of SN1998bw rather than SN1994I. All magnitudes have been 
corrected for foreground and host extinction. 
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Figure 9. GRB 090618: Rc - i for the afterglow, SN and host. 
The observed change in colour is not constant (as is expected for 
light from only a GRB) but increases with time. Such behaviour 
would be expected if a component of the flux was coming from 
a core-collapse SN. At late times the colour is that of the host 
galaxy {Rc - i ^ 0.22). 



brighter than SN1998bw in the V-band when host extinction 
is considered (see Table [7]). 

Additional evidence for an association of a SN with 
GRB 090618 is the change in Rc - i (of the afterglow, SN 
and host galaxy) over time (Figure [9|. At early times a blue 
colour is seen {Rc~i ~ 0.1), which we attribute to the after- 
glow. However, over time the colour index increases, which 
is not expected for an afterglow (i.e., solely synchrotron ra- 
diation), but is indicative of a component of hght coming 
from a core-collapse supernova. At late times as hght from 
the supernova fades away, the colour index decreases and 
approaches that of the host galaxy {Rc — i ~ 0.22, which 







Figure 10. GRB 090618: Rest-frame X-ray to optical SEDs: 
SEDl (t-to = 0.059 days); SED2 (t - to = 0.238 days); SEDS: 
{t-to = 0.568 days); SED4: {t - = 1.68 days). The results of 
the SED modelling are: (1) small rest-frame dust extinction, with 
Ay ~ 0.3 ± 0.1; (2) Each epoch is well fit by a broken power-law 
(i.e., with a cooling break between the optical and X-ray energy 
bands) with ^optical ~ 0.5 and feray = Poptical + 0.5; (3) the 
break frequency is clearly decreasing with time, which indicates 
an ISM environment rather than a wind environment. 



is typical of GRB host galaxies (e.g.. 
Table 8). 



Savaglio et al. 2009: 



3.2.3 The Spectral Energy Distribution 

The rest-frame X-ray to optical SEDs were determined for 
four epochs: t - to = 0.059, 0.238, 0.568 and 1.68 days. The 
times of their occurrence are highlighted on Figure [5] and 
are referred to in Figure [TO] as SEDl to SED4. For SEDl, 
where there is no contemporaneous X-ray and optical data, 
the X-ray data have been interpolated to the time of the op- 
tical observations. For all of the SEDs, the X-ray data have 
been derived from a broader interval so as to collect enough 
photons, where a normalization coefficient was determined 
and the interpolation performed through a power-law, which 
proved to be a good approximation to the data. 

The paucity of optical points prevents us from discrim- 
inating between the different optical extinction profiles (i.e., 
MW, LMC & SMC), thus we considered only the most suc- 
cessful template, i.e., an SMC-type template. The SEDs 
were fit using the SMC dust extinction profile from Pei 
(1992) and combined with a parametrisation of the pho- 
tometric cross-section to account for the soft X-ray absorp- 
tion. They were then applied to a range of models, in par- 
ticular a simple power-law and a broken power-law. The 
simple power-law case was immediately rejected. Taking the 
first epoch as an example (similar behaviour was also seen 
in the other three epochs), the slope of the optical points is 
Popt ~ 0.3 while that of the X-rays is Px ~ 1-0. These values 
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Table 6. Best-fitting results to SED modelling for GRB 090618 



SED 


T - To (days) 


/^optical 


0xray 


i^c (Hz) 


Av,reat (mag) 


xVdof 


SEDl 
SED2 
SED3 
SED4 


0.059 
0.238 
0.568 
1.680 


55+0 06 
0.48l»:« 
0.50l°;i? 


1 14+0.02 
^•1^^-0.01 

n QS+0 06 

1 oo+°i° 


3.9+li X 1017 
5.4+*'6 X 1016 

i.etit X 1016 

7.6l6'3 X 1015 


24+"-°^ 
n oQ+0.09 
"■•-'■J_0.09 
24+°'°^ 


1.0 
0.96 
0.92 
0.96 



GRB 090618 



1 1 1 




1 1 .J--- 1 . 
























SED4 - 






(1.68 days) 












Prudio 1 1 1 











le+08 le+09 k+10 k+11 le+12 le+13 le+14 le+15 le+16 le+17 le+18 le+!9 
v(Hi) 



Figure 11. GRB 090618: Rest-frame X-ray to radio SED at 
t — to = 1.68 days. The typical synchrotron frequency has been 
estimated from our data, where the optical and X-ray spectral in- 
dices have been found from the SED fitting to be: ^optical = 0-5; 
l^x = l3optical~^0-5- We have used the theoretically-expected value 
of Pradio = 1/3 and found Um = 3.66 X IQH Hz. The maximum 
flux was found to be Fmax = 1.43 mjy. 

are not consistent with a common value for p. If the optical 
points had a larger value of P than the X-ray, in principle 
the two values could be reconciled when considering dust 
effects. However the opposite case is seen here, indicating 
the only way to reconcile the two frequency regimes is by 
invoking a frequency break (i.e., cooling break, Vc) between 
the two. 

Thus the optical and X-ray spectral indices, as well as 
the local extinction were allowed to vary during the fit. Af- 
ter comparison of the initial results it was seen that the 
results of the fit for each epoch indicated a low value of 
the extinction, each with a slightly different value. We thus 
constrained the extinction within the same range for each 
epoch (as it is physically unlikely that the extinction would 
change from epoch to epoch), taking the best determined 
value from SED3, and allowing it to vary within its range 
(0.24 < Av < 0.42 mag) and re-fitted the other epochs. The 
best-fitting results are displayed in Table [S] 

The temporal evolution of the break frequency was also 
investigated (Table [6l Column 5). Sari, Piran & Halpern 
(1999) have shown that the cooling frequency evolves differ- 
ently depending on the geometry of the afterglow: fc oc t~^^^ 
for a spherical evolution, and ~ constant in a jet-like evo- 
lution. It was seen from the X-ray and optical light curves 
that a break occurred around t — to ~ 0.5 days, thus SEDl 
and SED2 occur before the break (i.e., spherical evolution), 
while SED3 and SED4 occur after the break (i.e., a jet-like 



evolution). For the first two epochs, a power-law was fit to 
the data (i.e., Vc = A t^). The minimum value of the 
power-law index that can be reasonably fit to the data is 
Vc oc i-0-72±o.i8 Idof = 2.08), which is consistent with 
that expected from theory. For the latter two epochs, the 
large uncertainties derived from the SED fitting (Table [6| 
do not allow us to investigate whether the frequency of the 
cooling break is still evolving or not, thus it can be constant 
as expected within the derived errors. 

The main conclusions we have drawn from the SED 
modelling are: (1) small rest-frame dust extinction, with 
Av ~ 0.3±0.1; (2) Each epoch is well fit by a broken power- 
law (i.e., with a cooling break between the optical and X-ray 
energy bands) with PopUcai ~ 0.5 and Pxray = PopUcai +0-5; 
(3) the break frequency is clearly decreasing with time, indi- 
cating an ISM environment rather than a wind environment 
(as one expects the break frequency to increase in a wind 
environment). The results of the SED fitting appear quite 
typical of general afterglow broadband spectra, as does its 
observed time evolution. 

3.2.4 The Jet-break, GRB energetics & Predicted Radio 
Light- Curves 

The rest-frame 1 — Iff* keV isotropic energy release in 
gamma-rays was determined by Ghirlanda, Nava & Ghis- 
ellini (2010) in their Table 1: E-^^i^o ~ 2.57 x lO'^^ ergs. 

A break was seen in the panchromatic LC (Figure[5]), oc- 
curring at both optical and X-ray wavelengths a,tt — to ~ 0.5 
days. If the break is interpreted in the context of a single jet 
(double sided with sharp edges), using a simple model based 
upon the standard theoretical framework (e.g., Rhoads 1997; 
Sari, Piran & Halpern 1999) and assuming: (1) forward- 
shock emission; (2) synchrotron radiation; (3) adiabatic evo- 
lution of the afterglow; (4) no self-absorption; (5) canonical 
medium density of n ~ 1 cm~'^; and (6) an isotropic equiva- 
lent kinetic energy comparable to that observed in gamma- 
rays, as well as the definition for Oj^t fi'om Sari, Piran & 
Halpern (1999), we estimate an opening angle of Oj^t ~ 1.5°. 
This in turn implies a corrected, observed gamma-ray emis- 
sion of E.y^g ~ 8.16 X 10''^ ergs. If the break is indeed due to 
a jet a homogeneous ISM is favoured, as a clear break is not 
expected for a wind environment (e.g., Kumar & Panaitescu 
2000). 

We have been fortunate with our contemporaneous ra- 
dio detection at 4.86 GHz at t — to = 1-68 days, and have 
modelled our radio, optical and X-ray data as the SED 
shown in Figure [TT] We have assumed a jet-like evolution of 
the afterglow for t — to > 0.5 days (i.e., f{t) oc t^i'^^), and 
neglected effects due to self- absorption, as well as assum- 
ing oc u^^^ for frequencies below the typical frequency 
(i.e., u < Um). From our SED modelling at t — to = 1.68 
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days we find Vsa < l^radio < Vm < l^opUcal < Vc < l^xray 

and have estimated the typical synchrotron frequency from 
our data to be Vm ~ 3.66 x 10^^ Hz and the peak flux to be 
Fmax = 1.43 mjy (Figure [TT|) . Using our above assumptions, 
as well as the derived energetics of the GRB, we have cal- 
culated the expected typical frequency of the synchrotron 
photons and compared it with the observationally-derived 
value. In the case where the magnetic field energy density 
is of order the electron energy density (i.e., eb ~ £e) we 
estimate eb ~ £e ~ 0.045. Alternatively, if eb ~ O.lSe, we 
estimate eb ~ 0.007 and Ee ~ 0.07. Additionally, if we use 
the value of the cooling/break frequency at t — to — 1.68 
days (I'c = 7.6tg'2 x 10^^ Hz) that was derived from the 
SED modelling to estimate eb and the medium density n, 
we derive eb ~ 0.001 if we assume n ~ 1. Alternatively, if we 
assume eb ~ 0.01, we derive a medium density n ~ 0.05 pro- 
ton cm~^. These micro-physical values are similar to the as- 
sumed values used in the literature (e.g., Zhang et al. 2006). 

Using our simple model we have neglected self- 
absorption effects. However, this assumption may not be 
completely reasonable, and the effect of self-absorption 
would result in a lower typical frequency and a higher peak 
flux. We calculate, using typical values of the micro-physical 
parameters as well as the definition for i/sa from Granot, Pi- 
ran & Sari (1999), a rest-frame self- absorption frequency of 
1/ « 2.25 X lO'^'' Hz, which is located in the radio frequency 
range of our detections. 

Assuming that the afterglow has a jet-like evolution af- 
ter ~ 0.5 days, we can estimate the electron energy distribu- 
tion index p from the optical and X-ray LCs by assuming t^^ 
(i.e., a = p) after the break (Sari, Piran & Halpern 1999). 
The optical data imply p « 1.6 and the X-ray data imply 
p ~ 1.75. These values are consistent with derived values of 
p before the break (i.e., for a spherical evolution of the fire- 
baU). Using the X-ray data before the break (i.e., t-3p/4+i/2 
above the cooling frequency), the electron index is found to 
be p « 1.72, while from the optical data before the break 
(i.e., t~3(p~i)/'* below the cooling frequency), the value of 
Q implies p « 1.87. While these values are lower than the 
predicted universal value of p ~ 2.2 — 2.3 (Kirk et al 2000; 
Achterberg et al. 2001), our derived value of p can be ac- 
commodated for when it is considered that p cannot have a 
single value for all GRBs (e.g., Kann et al. 2006; Shen et al. 
2006; Starling et al. 2008; Curran et al. 2010), but has a wide 
range of possible values (e.g., Curran et al. 2010 used X-ray 
spectral indices of GRB afterglows to parametrize the un- 
derlying distribution of p, finding the range consistent with 
a Gaussian distribution centered on p = 2.36 and having a 
width of 0.59. Our derived value of p ~ 1.7 — 1.8 falls at the 
lower end of this distribution). 

The initial Lorentz factor of the ejecta can be estimated 
by considering the deceleration time. Assuming that the on- 
set of the afterglow is masked by the steep decay component 
(i.e., the tail of the prompt emission), using the definition 
from Sari & Piran (1999), and correcting for observer time, 
we estimate an initial Lorentz factor of 70 > 100. 

The derived radio light curves are shown in Figure [51 
where the fiux scalings have been derived from the X-ray to 
radio SED at t — to — 1.68 days fFigure lTTj) . Both a spherical 
{fv{t) oc f^/^) and jet-hke (f„{t) cx t'^^'^) evolution of the 
afterglow are considered, and we conclude that all of the 
radio observations can be suitably explained by a jet-like 




Figure 12. Cumulative fraction plots for the two comparison sce- 
narios: {right) between all SNe and (left) considering only those 
events where an estimation of the host/rest-frame extinction has 
been made. In the latter case, the probabiUties that the three 
sets of SNe are drawn from the same parent population are (1) 
GRB/XRF SNe & All Ibc SNe: P = 0.88, (2) GRB/XRF SNe 
& only Ic SNe: P = 0.54, which further supports the idea that 
GRB-SNe and type Ic SNe have similar progenitors. 

evolution of the fireball, though we note that the scarcity of 
radio detections limits our ability to confidently determine 
the actual evolution of the afterglow beyond doubt. 

4 FURTHER DISCUSSION & CONCLUSIONS 

4.1 The SNe accompanying GRBs 060729 & 
090618 

We have presented photometric evidence for supernovae as- 
sociated with GRBs 060729 and 090618. For GRB 060729, 
peak apparent magnitudes of the associated supernova are 
Tie = 23.80 ± 0.08 and Jc = 23.20 ± 0.06. When rest-frame 
extinction was considered using the analysis performed by 
Schady et al. 2010 {E{B - F) 0.06 mag), the peak, 
rest-frame absolute V-band magnitude was shown to be 
Mv = —19.43 ±0.06 which is approximately the same peak 
brightness as SN1998bw in the V-band. 

For GRB 090618 the peak apparent magnitudes of the 
associated supernova are Rc = 23.45 ± 0.08 and Ic = 
23.00 ± 0.06. The rest-frame extinction local to the GRB 
was determined from the X-ray to optical SED fitting, and 
found to be Ay = 0.3 ± 0.1. In turn these values imply a 
peak, rest-frame absolute 1/-band magnitude of the SN as- 
sociated with GRB 090618 of Mv = -19.75 ± 0.13, which 
is 0.3 mag brighter than SN1998bw in the V-band. 

4.2 Comparison with existing GRB &: 
XRF-associated SNe 

To put our detections into context, we compiled from the 
literature peak, rest-frame 1^-band absolute magnitudes for 
two samples of type Ibc SNe, incorporating the values of 
the host/rest-frame extinction (taken at the location of the 
GRB or SN when available; i.e., Kann et al. 2006; Modjaz et 
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Table 7. Peak Rest-Frame V-band Absolute Magnitudes for GRB & XRF-producing SNe 



GRB 


SN 


Rcdshift (z) 


^V. f or e ground 


J^V,host " 


A/^""'' (mag)''-'= 


Reference 


GRB 970228 




0.695 


0.543 


0.15 


— 18.56 ± 0.30 


(1): (2), (3) 


GRB 980326 




s; 1 


0.26 




-19.5 


(4) 


GRB 980425 


1998bw 


0.0085 


0.18 


0.05 


— 19.42 ± 0.30 


(3), (5), (6), (7), (8), (31) 


GRB 990712 




0.434 


0.09 


1.67 


—20.22 ± 0.20 


(3), (10), (11), (12), (31) 


GRB 991208 




0.706 


0.05 


0.76 


— 19.46 ± 0.75 


(9), (16) 


GRB 0U0911 




1.058 


0.38 


0.20 


— 18.31 ± 0.15 


(9), (16) 


GRB 011121 


2001ke 


0.36 


1.33 


0.39 


— 19.59 ± 0.33 


(3), (13), (14), (16) 


GRB 020405 




0.698 


0.14 


0.15 


— 19.46 ± 0.25 


(3), (15), (16), (31) 


GRB 020410 




~ 0.5 


0.40 


0.0 


~ —17.6 


(3), (17) 


Arlr UzUyUo 




U.zoi 


n nn 

u.oy 


U.O 


1 Q Qo J- n Qn 
— io. oy ± U.oU 


(3), (18), (31) 


GRB 021211 


200211 


1.006 


0.08 


0.0 


-18.27 ±0.60 


(9), (19), (16) 


GRB 030329 


2003dli 


0.169 


0.07 


0.39 


-19.14 ±0.25 


(3), (16), (20), (31), (32) 


XRF 030723 




Si 0.4 


0.089 


0.23 


« -17.9 


(3), (9), (34) 


GRB 031203 


20031w 


0.1055 


2.77 


0.85 


-20.39 ±0.50 


(3), (21), (22), (31) 


GRB 040924 




0.859 


0.18 


0.16 


-17.47 ±0.48 


(23) 


GRB 041006 




0.716 


0.07 


0.11 


-19.57 ±0.30 


(3), (16), (24) 


GRB 050525A 


2005nc 


0.606 


0.25 


0.32 


-18.76 ±0.28 


(3), (25), (26), (33) 


XRF 060218 


2006aj 


0.033 


0.39 


0.13 


-18.76 ±0.20 


(3), (27), (28), (31) 


GRB 060729 




0.54 


0.11 


0.18 


-19.43 ±0.06 


This paper 


GRB 080319B 




0.931 


0.03 


0.05 


-19.12 ±0.40 


(29), (33) 


GRB 090618 




0.54 


0.27 


0.3 


-19.75 ±0.14 


This paper 


GRB 091127 


2009nz 


0.49 


0.12 


0.0 


-19.00 ±0.20 


(30) 



"Host extinction where available. 

Cosmological Parameters used: Ho = 71 km Mpc~^, f2jvf = 0.27, Ha = 0.73. 
'^Wherever errors are not quoted in the literature conservative errors of 0.4 mag are used. 
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al. 2008; Levesque et al. 2010), and applying the SMC red- 
dening law. The two samples are: (1) those associated with 
GRBs & XRFs, and (2) local type Ibc SNe. We note that 
we have limited this analysis to incorporate only GRB/XRF 
events where an optically-bright SN has been positively de- 
tected. These samples, along with their respective references, 
are listed in Tables [7] and [8| 

When comparing these two samples we are attempt- 
ing to answer the question: "Are the progenitors of GRB & 
XRF associated SNe the same as those of local type Ibc SNe 
without an accompanying GRB /XRF trigger?" by testing if 
the distribution of the peak magnitudes of the two sam- 
ples of supernovae are different. To do this we performed 
a Kolmogorov-Smirnov (KS) test on the two samples. First 
we compared the entire GRB/XRF sample (A'^ = 22) with 
all of the type Ibc SNe {N = 34), finding a modest proba- 
bility that the two samples are drawn from the same parent 
population of P = 0.12. When we compare the GRB/XRF 
sample with only the type Ic SNe {N — 19) we find a similar 
probability of P = 0.16. 

However, when we limited the samples to include only 
those SNe where an estimation of the host /rest-frame ex- 
tinction has been made, we see an increased probability 
between the datasets: P = 0.88 between the GRB/XRF 



SNe (iV = 21) and aU of the type Ibc SNe (iV = 17), and 
P = 0.54 between the GRB/XRF SNe and only the type 
Ic SNe (A^ = 12). For these samples a higher average peak 
magnitude is also seen among the local type Ibc SNe sam- 
ple, with the local type Ibc SNe having My, ibc = —18.93. 
This is an increase in average peak brightness of ~ 0.3 mag, 
and is likely due to the inclusion of five additional bright 
type lb events with known host-extinction (SN: 1983N, 
1983V, 1984L, 1985F & 1991D) with the local Ic SNe sam- 
ple (N=12). The average peak magnitude of these five SNe is 
Mv,ib = —19.34, which places them among the brightest of 
the type lb events. For comparison, a study by Richardson et 
al. (2002) analysed the absolute peak magnitudes of samples 
of all types of nearby SNe, and distinguished between "Nor- 
mal" and "Bright" type Ibc events, with the latter having 
Mvjbc ~ —19.72 ± 0.24, comparable with the average peak 
magnitudes of these five SNe. Thus the inclusion of these 5 
bright lb events has increased the average peak brightness 
of the type Ibc SNe sample, as well as increased the proba- 
bility of association between our samples of host extinction- 
corrected type Ibc SNe and the GRB/XRF-SNe. The results 
of the K-S test are summarized in Table |9] and cumulative 
fraction plots are shown in Figure [12] 

It is worth noting, however, that while the increase in 
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Table 8. Peak Rest-Frame V-band Absolute Magnitudes for Local type Ibc &l Ic SNe 


Type 
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Redshift {z) 
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M^^""^ (mag)'''^ 
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1.18 


— 1 Q 1 2 + 41 




ID 




U.UlU / 


'i'i 
V.oo 




17 Kfi -1- 40 


(Q\ 


lb 


1984L 


U.UUOZOl 


OS 
u.uo 




u.u 


1 S R4 -1- 40 


(^ 0"! 


TK 




n nni ^7 


Of, 
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Ic 
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Ic 
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007'^l S 
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2.53 
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lb 
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Ic 


1991N 
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1999ex 
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- 


-17.67 ±0.26 


(23) 


lb 
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0.39 




-17.13 ±0.40 


(24) 


Ic BL 


2002ap 
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0.29 


0.0 


-17.73 ±0.21 


(2), (25) 


Ic 


2003L 
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-18.90 ±0.40 
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0.018826 


0.14 


0.29 
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"Host extinction where available. 

''Cosmological Parameters used: Ho = 71 km s~^ Mpc~^, Q,m = 0.27, Ha = 0.73. 
'^Wherever errors are not quoted in the literature conservative errors of 0.4 mag are used. 
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Table 9. Kolmogorov-Smirnov test results 



Datasct 


Number of Data points 


Mean 


Standard Deviation 


pa 




comments 


GRB/XRF-associated SNe 


22 


-19.02 


0.77 






all events 


Local type Ic SNe 


19 


-18.73 


1.00 


0.16 


0.33 


all events 


Local type Ibc SNe 


34 


-18.59 


1.04 


0.12 


0.31 


all events 


GRB/XRF-associated SNe 


21 


-19.00 


0.78 






only those with host Ay 


Local type Ic SNe 


12 


-18.75 


1.03 


0.54 


0.27 


only those with host Ay 


Local type Ibc SNe 


17 


-18.93 


0.97 


0.88 


0.18 


only those with host Ay 



"Probability and maximum difference between the GRB/XRF SNe sample and the local SNe sample. 
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probability may really be due to an increased association 
between the datasets, the caveat of smaller sample sizes is 
that as one goes to smaller samples it is harder to obtain a 
statistically significant discrepancy. 

Our analysis is not unique, with previous studies hav- 
ing performed similar analysis, which our results generally 
support. Richardson 2009 (R09) found from a sample of 14 
GRB/XRF-SNe with (mostly) known values of the host ex- 
tinction an average Mv,peak = —19.2 ± 0.2, with a = 0.7, 
which our results are in agreement with. R09 also com- 
pared the GRB/XRF-SNe sample with a sample of stripped- 
envelope SNe, which included types lb, Ic and lib finding the 
GRB/XRF-SNe sample brighter by ~ 0.8 mag. We find that 
our complete GRB/XRF-SNe sample is ~ 0.4 mag brighter 
than the complete local type Ibc SNe (but only 0.1 mag 
when only events with known host extinction are consid- 
ered), which is somewhat less than R09, however our sample 
does not include any type lib SNe. 

Ferrero et al. 2006 (F06) undertook a slightly different 
study, comparing the existing GRB/XRF-SNe sample (to 
date) with that of XRF 060218. They furthered an original 
analysis performed by Zeh et al. (2004), and calculated for 
each GRB-SN event the luminosity ratio k and stretch fac- 
tor s in comparison with SN1998bw. They find for their host 
extinction-corrected sample of GRB/XRF-SNe clustering in 
the range 0.6 < k < 1.5, implying a range of absolute mag- 
nitudes at peak brightness of -1-0.55 mag < My^"" < -0.44 
mag, which is in good agreement with our results. They 
conclude that the width of the GRB/XRF-SN luminosity 
function is at least 2 mag, and there was no evidence that 
the luminosity function evolved with redshift. Our analysis 
concurs with both of these results. 

In conclusion, while the complete sample of GRB/XRF 
SNe are generally brighter (Mv,grb = -19.02, a = 0.77) 
than the complete local type Ibc SNe {Mvjtc ~ —18.59, 
a = 1.04) and the type Ic SNe sample {Myjc = -18.73, 
cr = 1.00), our test does not address factors such as host and 
progenitor metallicity and typical outfiow velocities. Thus 
we cannot rule out the null-hypothesis that the samples of 
SNe are drawn from the same parent population. 
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